THE LEVEL SET METHOD FOR THE TWO-SIDED MAX-PLUS 

EIGENPROBLEM 



STEPHANE GAUBERT AND SERGEI SERGEEV 

Abstract. We consider the max-plus analogue of the eigenproblem for matrix pencils, 
A ® x = A ® B (g> x. We show that the spectrum of (A, B) (i.e., the set of possible values 
of A), which is a finite union of intervals, can be computed in pseudo-polynomial number 
of operations, by a (pseudo-polynomial) number of calls to an oracle that computes the 
value of a mean payoff game. The proof relies on the introduction of a spectral function, 
which we interpret in terms of the least Chebyshev distance between A®x and \ ®B®x. 
The spectrum is obtained as the zero level set of this function. 



1. Introduction 

1.1. Motivations and general information. Max-plus algebra is the analogue of linear 
algebra developed over the max-plus semiring which is the set R max = MU{— 00} equipped 
with the operations of "addition" a@b := aV b = max(a, b) and "multiplication" a®b := 
a + b. The zero of this semiring is —00, and the unit of this semiring is 0. Note that 
a -1 in max-plus is the same as —a in the conventional notation. The operations of the 
semiring are extended to matrices and vectors over M max . That is if A = (a y -), B = (6^) 
and C = ( matrices of compatible sizes with entries from M max , we write C = A\/ B 

if Cij = a,ij V bij for all i,j and C — A <g) B if Cy = v k{ a ik + bj~j) for all 

We investigate the two-sided eigenproblem in max-plus algebra: for two matrices A,B& 
M™, find scalars A 6 M ma x called eigenvalues and vectors x 6 M™^ called eigenvectors, 
with at least one component not equal to —00, such that 

(1) A®x = \®B®x, 
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where the operations have max-plus algebraic sense. In the conventional notation this 
reads 

(2) max(ajj + Xj) = A + max(6jj + xA for i = 1, . . . , m. 

i=i ' i=i 

The set of eigenvalues will be called the spectrum of (A, B) and denoted by spec(A, B). 

When B is the max-plus identity matrix / (all diagonal entries equal and all off- 
diagonal entries equal — oo), problem ([1]) is the max-plus spectral problem. The latter 
spectral problem, as well as its continuous extension for max-plus linear operators, is of 
fundamental importance for a wide class of problems in discrete event systems theory, 
dynamic programming, optimal control and mathematical physics [30], l3Tj . 

Problem ([T]) is related to the Perron- Frobenius theory for the two-sided eigenproblem in 
the conventional linear algebra, as studied in [3lJ[35]. When both matrices are nonnegative 
and depend on a large parameter, it can be shown following the lines of [U Theorem 1] 
that the asymptotics of an eigenvalue with nonnegative eigenvector is controlled by an 
eigenvalue of (JTJ. This argument calls for the development of two-sided analogue of the 
tropical eigenvalue perturbation theory presented in [3J 12] • 

A specific motivation to study the two-sided max-plus eigenproblem arises from discrete 
event systems. In particular, systems of the form A®x = B®a;orA®i^B®i appear 
in scheduling. Indeed, when A = 0, the system of constraints (EJ can be interpreted in 
terms of rendez-vous. Here, Xj represents the starting time of a task j (for instance, 
the availability of a part in a manufacturing system). The expression max" =1 (ajj + Xj) 
represents the earliest completion time of a task which needs at least time units to 
be completed after task j started. Thus, the system A <g> x = B <g> x requires to find 
starting times such that two different sets of tasks are completed at the earliest exactly 
at the same times. In many situations, such systems cannot be solved exactly, and a 
natural idea is to calculate the minimal Chebyshev distance between A £g> x and B (g) x. 
Theorem 4 below determines this minimal distance. It may be also of interest to solve 
perturbed problems like A®x = \®B®x, as in ([T]). Such problems express no wait 
type constraints. Indeed, y := A <g> x and z := B £g> x may be thought of as the outputs 
of two different systems A and B, with a common input x. The time offsets between 
output events are represented by the differences y^ — yj, for all i,j, where the difference is 
understood in the usual algebra (these quantities belong to the "second order" max-plus 
theory, see e.g. [IT]). No wait constraint may require that — yj take prescribed values, 
for some pair The condition that y = A <g) z, i.e., y = A + z, for some A G R, 

means precisely that the time offsets are the same for the two outputs y and z. Hence, an 
input x solving A®x = A®-B®x has the property of making A and B indistinguishable 
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from the point of view of no wait output constraints. An example of such a situation is 
demonstrated on Figure [2] in Subsection 13.21 

Problems of a related nature, regarding the time separation between events, arose for 
instance in the work of Burns, Hulgaard, Amon, and Borriello [12], following the work 
of Burns on the checking of asynchronous digital circuits [13] . Moreover, systems of the 
form A <g> x ^ B <g> x represent scheduling problems with both AND and OR precedence 
constraints, studied by Mohring, Skutella, and Stork [57] . 

Similar motivations led to the study of min-max functions by Olsder [39J and Gunawar- 
dena [29J . Such functions can be written as finite infima of max-plus linear maps, or finite 
suprema of min-plus linear maps. They also arise as dynamic programming operators of 
zero-sum deterministic games. In particular, the fixed points and invariant halflines of 
min-max functions studied in [TBI [25] can be also used to compute values of zero-sum 
deterministic games with mean payoff [25J H3] . A correspondence between the computa- 
tion of the value of mean payoff games and two-sided linear systems in max-plus algebra 
has been established in [5]; we shall exploit here the same correspondence, although in 
different guises. 

In max-plus algebra, a special form of min-max functions appears in Cuninghame- 
Green [19], under the name of AA*-products. The same functions appear as nonlinear 
projectors on max-plus cones playing essential role in the max-plus analogue of Hahn- 
Banach theorem [18]. [33] . The compositions of nonlinear projectors are more general min- 
max functions, and they appear when one approaches two-sided systems A £g> x = B <g> y 
and A <g> x = B <g> x |20j , and intersections of max-plus cones [281 H2] • It is immediate to 
see that (pQ) is a parametric version oiA®x = B®x. 

In max-plus algebra, partial results for Problem ([T]) have been obtained by Binding and 
Volkmer [10] , and Cuninghame- Green and Butkovic [151 121] • m particular, Cuninghame- 
Green and Butkovic [2TJ give an interval bound on the spectrum of (|T}) in the case 
where the entries of both matrices are real. Besides that, both papers treat interesting 
special cases, for instance when A and B square, or one of them is a multiple of the other. 

The spectrum of ([T]) is generally a collection of intervals on the real line. By means 
of projection, this follows from a result of De Schutter and De Moor [ID] that solution 
set to the system of max-plus (in) equalities is a union of convex polyhedra. Note that 
the approach of |40j, related to Develin-Sturmfels cellular decomposition |22j, can be also 
used for solving A®x = A £g> B®x and more general problems of max-plus linear algebra. 

1.2. Contents of the paper. In the present paper, we first show that (fT]) can be viewed 
as a fixed-point problem for a family of parametric min-max functions h\. Based on this 
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observation, we introduce a spectral function s(A) of flTJ, defined as the spectral radius of 
h\. The zero level set of s(A) is precisely spec(A, B). More generally, s(A) has a natural 
geometric sense, being equal to the inverse of the least Chebyshev distance between A <g> x 
and A <8> B <g> x. 

The function s(A) is piecewise-affine and Lipschitz continuous, and it has an affine 
asymptotics at large and small A. In an important special case when none of the matrices 
A and B have an identically — oo column, the asymptotics is just A + «i at small A, and 
— A + ct2 at large A, in the conventional arithmetics. We also give bounds on the spectrum 
of two-sided eigenproblem, which improve and generalize the bound of Cuninghame-Green 
and Butkovic [T51 [2TJ. In the case when the entries of A and B are integer or — oo, this 
allows us to show that all affine pieces of s(A) can be identified in a pseudopolynomial 
number of calls to an oracle which identifies s(A) at a given point. Importantly, s(A) 
can be interpreted as the greatest value of the associated parametric mean-payoff game 
and it can be computed by the policy iteration algorithm of [T6| [23] . as well as by the 
value iteration of Zwick and Paterson [H] or the subexponential method of Bjorklund 
and Vorobyov [IT]. This leads to a procedure for computing the whole spectrum of p]). 
To our knowledge, no such general algorithm for computing the whole spectrum of ([I]) 
was known previously. We also believe that the level set method used here, relying on 
the introduction of the spectral function, is of independent interest and may have other 
applications. 

In some cases the spectral function can be computed analytically. In particular, we 
will consider an example of [41J, where it is shown that any finite system of intervals and 
points on the real line can be represented as the spectrum of ([T|). 

The paper is organized as follows. In the remaining subsection of Introduction we 
explain the notation used in the rest of the paper. In Section 2 we consider two-sided 
systems A ® x = B £g> y and A <g> x = B <g> x. We relate the systems A ® x = B ® x to 
certain min-max functions and show that the spectral radii of these functions are equal 
to the inverse of the least Chebyshev distance between A <g> x and B®x. In Section 3, we 
introduce the spectral function of two-sided eigenproblem as the spectral radius of a nat- 
ural parametric extension of the min-max functions studied in Section 2. We give bounds 
on the spectrum of two-sided eigenproblem and investigate the asymptotics of s(A). We 
reconstruct the spectral function and hence the whole spectrum in a pseudopolynomial 
number of calls to the mean-payoff game oracle. 

1.3. Notation. For the sake of simplicity, the sign £g> will be usually omitted in the 
remaining part of the paper, or even replaced with + if scalars are involved. In particular 
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we write Ax for A <g) x and A + x for A (g> x, where A G I^max™ ( a matrix), x G K max (a 
vector) and A G K m ax (a scalar). Moreover in the remaining part of the paper we will 
prefer conventional arithmetic notation: the four arithmetic operations a + b, a — b, ab 
and a/6 on the set of real numbers (scalars) will have their usual meaning. However we 
often use V for max and A for min (also componentwise). The actions of max-plus linear 
operators, their min-plus linear residuations and nonlinear projectors onto max-plus cones 
(defined below in Subsection 12. ip . which will appear as Ax, A^y, Pxy etc., should not be 
confused with any conventional linear operator. The notations like A^B or PaPb should 
be understood as compositions of the corresponding operators rather than any kind of 
matrix multiplication between them (in the case of A*B above, B is max-plus linear and 
v4" is min-plus linear). 

Such notation is implied by the methodology of the paper: we consider a problem of 
max-plus algebra, using nonlinear Perron-Frobenius theory and elementary analysis of 
piecewise-afline functions, where the max-plus arithmetic notation is not required or is 
inconvenient. The max-plus matrix product notation g), especially when mixed with the 
dual min-plus £g/, is also inconvenient when several compositions of such operators appear 
in the same formula or equation. 

The notation used in our main subsections is not new. Being different from Baccelli et 
al. [9] or Cuninghame-Green [19] , it basically follows Gaubert and Gunawardena [25 . 



2. Two-sided systems and min-max functions 

2.1. Max-plus linear systems and nonlinear projectors. Consider the m-fold Carte- 
sian product Mmax equipped with operations of taking supremum u V v and scalar "multi- 
plication" (i.e., addition) X®v := X + v. This structure is an example of semimodule over 
the semiring M max defined in the introduction. The subsets of 1™^ closed under these 
two operations are its subsemimodules. We will call them max-plus cones or just cones, 
by abuse of language. Indeed, there are important analogies and links between max-plus 
cones and convex cones [T8| |22| [26| W2\ . We also need the operation of taking infimum 
which we denote by inf or A. 

With a max-plus cone V C M^ ax we can associate an operator P v defined by its action 

(3) P v z = V{y G V | y ^ z}. 

Consider the case where V C M™ ax is generated by a set S G 1™^, which means that it is 
the set of bounded max-plus linear combinations 

(4) v = V X y + y. 
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In this case 

Pyz = V zjy + y, where 

(8) 

zjy = max{7 | 7 + y ^ z} = A (zj - = A (zj - 

jGsupp(j/) J = l 

with the convention (—00) + (+00) = +00. Here and in the sequel supp(y) := {i \ yi 7^ 
—00} denotes the support of y. Note that zjy = 00 if and only if y = — 00. 

Further we are interested only in the case when V is finitely generated. Let S = 
{y 1 , . . . , y n }, and 7] denote the set of indices where the minimum in zjy 1 is attained. The 
following result is classical. 

Proposition 1 (0, [TJJ ED])- Let a cone V C M™ ax fre generated by y 1 , . . . ,y n and let 

z G I^max- ^ e following statements are equivalent. 

1. zeV. 

2. P v z = z. 

3. LX=i T * = suppz. 

We note that the set covering condition 3. has been generalized to the case of Galois 
connections [6]. 

By this proposition, the operator Py is a projector onto V. It is an isotonic and +- 
homogeneous operator, meaning that z 1 ^ z 2 implies Pyz 1 ^ Pyz 2 , and that Py(\ + z) = 
A + Pyz. However, in general it is neither V- nor A-linear. 

A finitely generated cone can be described as a max-plus column span of a matrix 
A e R™™: 

n 

(6) span(A) := { V X { + A. { | X { G M max , i = l,...,n}. 

i=l 

In this case we denote Pa '■= -P S pan(A), and there is an explicit expression for this operator 
which we recall below. 

We denote R max := K m ax U {+00} and view A G M max as an operator from M max to 
M max . The residuated operator A^ from to M.™ is defined by 



rn 



(7) (A*y)j = yjA. 3 = A {-a t] + Vi ), 

i=i 

with the convention (— 00) + (+00) = +00. Note that this operator, also known as 
Cuninghame- Green inverse, sends M max to M™ ax whenever A does not have columns equal 
to —00. The term "residuated" refers to the property 

(8) Ax < y <^> x < A^y, 
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where ^ is the partial order on ]R™ ax or M max . Using (jSJ) we obtain 



(9) P A (z) = y(zjA.i)+A4 = AAl 



z. 



i=l 



In this form fl9]), the nonlinear projectors were studied by Cuninghame-Green [19] (as 
AA*-products). 

Finitely generated cones are closed in the topology induced by the metric 



(10) d(x, y) = max \e Xl — e 



in i 



1 1 



which coincides with Birkhoff's order topology. It is known [18, Theorem 3.11] that the 
projectors onto such cones are continuous. 

The intersection of two finitely generated cones can be expressed in terms of two-sided 
max-plus linear systems with separated variables Ax = By, by the following proposition. 



Proposition 2. Let A e M™^ 1 and B e K™ a x x n2 - 

1. If(x, y) satisfies Ax = By ^ — oo then z = Ax = By belongs to span(74)flspan(I?). 
Equivalently, PaPbz = PbPaz = z. 

2. If P a Pbz = — oo then there exist x and y such that Ax = By = z. 

This approach to two-sided systems is also useful in the case of systems with non- 
separated variables Ax = Bx, which is of greater importance for us here. This system is 
equivalent to 

Cx = Dy, where 

-(;)■ 

and I m = (5ij) G R^ax™ denotes the max-plus m x m identity matrix with entries 

0, if % = j, 
-oo, if i ^ j. 

In this case we have the following version of Proposition |2j 

Proposition 3. Let A, B e M™ ax n . 

1. If x satisfies Ax = Bx ^ — oo, then v = (z z) T , where z = Ax = Bx, belongs to 
span(C) fl span(.D). Equivalently, PqPdv = PdPc v — Pc v — v - 

2. If v = (z z) T 7^ — oo and Pqv = v, then there exist x such that Ax = Bx = v. 



(12) 8. 
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Pairs (x, y) 7^ — oo such that Ax = By = — oo are described by: Xi 7^ — oo A, = —00 
and yj 7^ —00 B.j = —00. Analogously, vectors x 7^ —00 such that Ax = Bx = —00 
are described by X{ 7^ —00 <^ A* = B.i = —00. Any such pair of vectors can be added 
to any other pair (x',y') or, respectively, vector x', and the resulting pair of vectors will 
satisfy the system if and only if so does (x', y') or, respectively, x'. Therefore, we can 
assume in the sequel without loss of generality that there are no such solutions, i.e., that 
1) A and B do not have —00 columns in the case of separated variables, 2) A and B do 
not have common —00 columns in the case of non-separated variables. 

2.2. Projectors and Perron-Frobenius theory. Suppose that a function / : R max — > 
R max is homogeneous, isotone and continuous in the topology induced by (fTUl) . As x i-» 
exp(x) yields a homeomorphism with R™ endowed with the usual Euclidean topology, we 
can use the spectral theory for homogeneous, isotone and continuous functions in R™ . We 
will use the following important identities, which follow from the results of Nussbaum |38j . 
see [5j Lemma 2.8] for the proof. 

Theorem 1 (Coro. of [35], Lemma 2.8]). Let f denote an order-preserving, additively 
homogeneous and continuous map from (R U {— oo}) n to itself. Then it has a largest 
eigenvalue 

r(f) := max{A | 3a; G R max , x ^ -00, A + x = f(x)}, 
which coincides with 

(13) r(/) = max{A | 3x e R max , x £ -00, A + x ^ f{x)}, 

(14) r(f) = inf{A | 3x G R n , A + x ^ f(x)}. 

Note that (Tl4l) is nonlinear generalization of the classical Collatz-Wielandt formula 
[55] . Equations f[T3T) and f[T4|) are useful in max-plus algebra, since they work for max- 
plus matrix multiplication as well as for compositions of nonlinear projectors. For (Till) it 
is essential that it is taken over vectors with real entries, and that the infimum may not 
be reached. Using (fT41) we obtain that the spectral radius of such functions is isotone: 
f(x) ^ g{x) for all x G R n implies r(f) ^ r(g). We next recall an application of (fl^|) to 
the metric properties of compositions of projectors, which appeared in [2H]- The Hilbert 
distance between u,v G R max such that supp(n) = supp(t>) is defined by 

(15) du(u, v) = max (u» — v 4 + Vj — Uj). 

jj'Ssupp(t)) 

If span(u) 7^ span(i>) then we set dn(u,v) = +00. Using fTl5l) we define the Hilbert 
distance between cones span(v4) and span(5), for A G R™ ax ni and B G E™^ 12 : 

(16) du(A,B) := min{<iH(^,^) I u G span(y4), v G span(i?), supp(w) = supp(w)}. 
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Theorem 2 (cp. [28], Theorem 25). Let A G E: 




(17) 



t{P A Pb) 



r(P B P A ) 



If dft(A,B) is finite then it is attained by any eigenvector u of PaPb with eigenvalue 
t(PaPb), and its image v by Pb- 

Proof. As supp(Pa-Pbm) ^ supp(Pew) C supp(w), it follows that PaPb and also PbPa 
may have finite eigenvalue only if span(A) and span(5) have vectors with common sup- 
port. This shows the claim for the case dfi(A, B) = +oo. 

Now let du(A,B) be finite. We show that — dn(u, v) = —dn(A,B) = t(PaPb)- Take 
arbitrary vectors u G span(A) and v G span(£?) with supp(w) = supp(t>), and let P u , resp. 
P v , be projectors onto the rays U = {A + u, A G M m ax}, resp. V = {A + v , A G M m ax}- 
As U C span(yl) and V C span(5), we have that P u ^ Pa and P v ^ Pg, hence P U P V ^ 
PaPb and, by the monotonicity of the spectral radius, r(P u P v ) ^ t{PaPb)- It can be 
shown that — f) is the only finite eigenvalue of P U P V , hence —dn(u,v) = r(P u P v ), 
and consequently —du(u,v) ^ t-^PaPb) and — cfn(A -6) ^ ^(PaPb)- Now observe that 
— c?h(w> v) = r(PuPv) is equal to the eigenvalue r{P A PB)- This completes the proof. □ 



In the case of the systems with non-separated variables, we will be more interested in 
the Chebyshev distance. For u, v G 1™^ with supp(u) = supp(f) it is defined by 



There is an important special case when Hilbert and Chebyshev distances coincide. 

Lemma 1. Let u, v G R™ ax be such that u ^ v and U{ = V{ for some i G {1, . . . , n}. Then 
dB_{u,v) = doo(u,v). 

Proof. First note that both du(u,v) and doo(u,v) are finite if and only if supp(w) = 
supp(f ). Only this case has to be considered. 

If u ^ v then \uj — Vj + Vi — u{\ ^ max(u-,- — Vj,ui — v{) for any j and I, hence 
d H (u, v) ^ d^u.v). 

Fixing I = i (assuming that Ui = v{) we obtain Uj — Vj + vi — ui = Uj — Vj and fixing 
j = i we obtain Uj — Vj +vi — ui — V\ — u\. Taking maximum over such terms only yields 
doo(u, v), hence d^u, v) ^ rf H (u, v). □ 



□ 



(18) 



max | Ui — Vi | . 

i£supp(i;) 



□ 
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Theorem 3. Let A, B G R££T, and let C and D be defined as in ([II]). Then 

(19) r{P c P D ) = r{P D P c ) = - min d^Ax.Bx). 

xtii^ max 

Proof. Theorem [2] implies that 

(20) r(PcPo) = — min{d}i(u, v) \ u G span(C),t> G span(.D).} 

Let u G span(C) and denote by P u the projector onto U := {A + w | A G M ma x}- Then u 
is an eigenvector of P u Pd which corresponds to the spectral radius of this operator, and 
applying Theorem [2] to the max cones U and span(D) we see that 

(21) d^u^Pou) = mm{du(u,v) \ v G span(D)}. 

Note that (1211) also holds if there is no v G span(D) with supp(u) = supp(f), in which 
case dn(u, Pdu) = +oo. This implies 

(22) r(P c P D ) = - min{d H (w, P D u) \ u G span(C)}. 
Observe that 

( Ax\ (AxNBx 
^ ' U ~ \Bx J ' ° U ~ \Ax A Bx 

for some x G R™ ax , and also that u and Pdu satisfy the conditions of Lemma [D unless 
Pdu = — oo. Hence dn(u,Prju) = d^u, Pdu) = d^Ax, Bx). Conversely, d^Ax, Bx) 
equals d-^(u, Pdu) for u = (Ax Bx) T . Hence the r.h.s. of (fT9|) is the same as the r.h.s. of 
(|22l) . which completes the proof. □ 

□ 



mx«2 

max " ^- ""max 



2.3. Min-max functions and Chebyshev distance. Let A G M^ax" 1 an< ^ B e 
In order to find a point in the intersection of span(A) and span(-B) (or equivalently, solve 
Ax = By), one can compute the action of (PaPb) 1 , for Z = 1, 2, . . . , on a vector z G M™ ax . 
Dually one can start with a vector a; G IR m a X and compute 

(24) x k = A^BB^Ax^ 1 , k ^ 1. 

We can assume that A and B do not have columns equal to — oo so that A*z G l™ 1 ^ and 
Bh G M^ ax for any z G M™ ax . 

If at some stage x k = x k ~ 1 ^ — oo then we can stop, x k is a solution of the system. If all 
coordinates of x k are less than those of x° then we can stop, the system has no solution. 
More details on this simple algorithm called alternating method can be found in [20] and 
[12], see also [7]. In particular, it converges to a solution with all components finite in a 
finite number of steps, if such a solution exists. 
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Let A,B G l'™^ 1 . A system Ax = Bx can be written equivalently as Cx = Dy with C 
and D as in ffTT]) . Applying alternating method (T24"|) to this system, i.e., substituting C 
and D for A and I? in ( 124)) we obtain x fc = g(x k ~ r ), where 



(25) g(x) = A^Ax A B^Bx A A*Bx A P fl A 



,r. 



As it is assumed that A and B do not have common — oo columns and hence C (and D) 
do not have — oo columns, g(x) G M™ ax for all x G fl^ax- 
It can be shown that (see also 



(26) r(g) = <^> Ax = Bx is solvable. 

In particular, if x is a fixed point of g then it satisfies Ax = Bx. For the function 

(27) /(x) = jA A^Bx A £ 8 Ax 

which appears in [23], it is also true the other way around, since 

Ax = Bx <^> Ax ^ Bx k Bx ^ Ax 

(28) 5Mi ^ x & A^Px ^ x <^> 

x A A^Bx A P tt Ax = x. 

We also introduce the function h: 

(29) h(x) := A tt Px A B % Ax. 

Although /, g and h are different functions, they have the same spectral radius, equal 
to the inverse minimal Chebyshev distance between Ax and Bx. To show this, we use 
the following identity. 

(30) — d^iu, v) = max{A : X+u^vkX+v^ u}. 

Theorem 4. Let A, B G R^ n - For C,D defined by (TTTT) . and f, g and h defined by 

dZZD, (ESD and (J235, 

(31) r(P c P D ) = r(P D P c ) = r(/) = r{g) = r(h) = - mm d^Ax, Bx). 

Proof. If f is an eigenvector of PdPc with a finite eigenvalue, then 0v is an eigenvector of 
g and P^v is an eigenvector PqPd, both with the same eigenvalue. The other way around, 
if x is an eigenvector of g with a finite eigenvalue, then (Ax Bx) T is an eigenvector of PdPc 
with the same eigenvalue. This argument shows that 1) either the spectral radii of P^Pc, 
PcPd and g are all finite or they all equal — oo, 2) the equality r(g) = r(PoPc) = ^(PcPd) 
holds true both in finite and in infinite case. 
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We show the remaining equalities. By (1131) . r(h) is the maximum of A which satisfy 

(32) 3a; e M max : A + x < A s 5a; A 
This is equivalent to 

(33) 3a; G R™ ax : A + Ar < 5a; & A + Bx < Aa; 
Using ( 1301) we obtain 

(34) r(h) — max — (^(Ac, £>a;) = — min ^(Aa;, 5x). 

It follows in particular that r(h) ^ and moreover, A ^ for any x satisfying (133]) . 
Applying (fl3|) to / and g we obtain that both r(f) and r(g) are equal to the maximum 
of A which satisfy 

(35) 3a; G M™ ax : A^O & X + Ax^Bx & A + < Ac 

As the first inequality follows from the other two, we obtain r(f) = r(g) = r(h). □ 

□ 

Functions /, g and h as well as projectors onto finitely generated max-plus cones and 
their compositions, belong to the class of min-max functions.Such functions were originally 
considered by Olsder [31] and Gunawardena [29] . See [TB] for a formal definition. In a 
nutshell, these are additively homogeneous and order preserving maps, every coordinate 
of which can be represented as a minimum of a finite number of max-plus linear forms, 
or as a maximum of a finite number of min-plus linear forms. It is important that any 
min-max function q : M" ax — > K max can be represented as infimum of finite number of 
max-plus linear maps Q^' meaning that 

(36) q(x) = AQ (p) x, 

v 

in such a way that the following selection 'property is satisfied: 

(37) Vx Bp : q(x) = Q {p) x. 

Note that taking infimum or supremum of vectors does not necessarily select one of them, 
and that selection property ([57]) is useful, e.g., for the policy iteration algorithm of [23J). 

In connection with the mean payoff games [23j [5], each matrix Q^' corresponds to a 
one player game, where the player Min has chosen her strategy and the player Max is 
trying to win what he can. 
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In particular, f(x),g(x) and h(x), respectively, are represented as infima of the max- 
plus linear maps F&\ and whose rows are taken from the max-plus linear forms 
appearing in ([27]) . ( 125]) and (129]) . respectively, in the following way: 



(38) F 2 



(p) 



4, 



—a ki + B k ., G i 
-b ki + A k .. 



(p) 



-Q>m + A k . } 
-bki + B k ., 
-Uki + B k ., 
-bki + A^. 



H (P) 




Here J,, denotes the ith row of the max-plus identity matrix, and the brackets mean that 
any possibility, for any k = 1, . . . , m and a k i ^ — oo or bki ^ — oo, can be taken (assumed 
that A and B do not have common — oo columns). Applying Collatz-Wielandt formula 
( I14p we obtain the following, some variants of which appeared in several contexts. 



Proposition 4 (Compare with [T6] 125] 124"] [8]). Suppose that a min-max function q : 
^max — Cm ^ s represented as infimum of max-plus linear maps 
selection property is satisfied. Then 

(39) r(g) = minr(Q (0 ). 

Proof. The spectral radius is isotone, hence r(q) ^ r(Q"') for all /. Using ( IT41 we conclude 
that for any e there is x G M m such that q(x) ^ r(g) + e + x. As g(x) = Q^x for some / 
and there is only finite number of matrices Q"\ there exists I such that 

(40) r(q) = inf{/i | 3x G R n , Q {l) x ^ \i + x} = r(Q w ). 
The proof is complete. □ 

□ 

Proposition H] can be derived alternatively from the duality theorem in (23 Theorem 19] 
(see also |24J). It is related to the existence of the value of stochastic games with perfect 
information |32j. Indeed, the spectral radius can be seen to coincide with the value of a 
game in which Player Max chooses the initial state, see [5] for more information. 

The greatest eigenvalue r{Q^>) of the max-plus matrix = (q$) G can be 

computed explicitly. It is equal to the maximum cycle mean of defined by 

(0 i (0 i (0 

a- ■ + a- ■ +... + (/•■ 

/ a 1 \ "1*2 "2*3 "fr*l 

(41) max max ^- — . 

l^k^n K 

This result is fundamental in max-plus algebra, see jU El [151 SO] for more details. 
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3. The spectrum and the spectral function 

3.1. Construction of the spectral function. Given A £ M™ ax n and B £ we 

consider the two-sided eigenproblem which consists in finding eigenvalues A £ M max and 
eigenvectors x £ M max (which have at least one component not equal to — oo), such that 

(42) Ax = A + Bx. 

The set of eigenvalues is called the spectrum of (A, B) and denoted by spec(A, B). 

Below we assume that A and B do not have — oo rows and common — oo columns. Note 
that the assumption about — oo rows can be made without loss of generality when the 
solvability of f|42|) is considered. Indeed, if the ith row of B is — oo then all variables Xj 
such that ^ — oo must be equal to — oo. Eliminating these variables as well as the 
corresponding columns in A and B and the ith equation, we obtain a new system where 
A or B may have — oo rows. Proceeding this way we either cancel the whole system in 
which case it is unsolvable, or we are left with a system where A and B (what remains of 
them) do not have — oo rows. This procedure can be run in 0(m 2 n) operations. 

The case of A = — oo appears if and only if A has — oo columns, and the corresponding 
eigenvectors are described by Xi ^ — oo <^ A4 = —00. In the sequel we assume that A is 
finite. 

Problem (p]) is equivalent to C(X)x = Dy, where C(A) £ and D £ M^ m are 

defined by 





( 43 ) C W = L,n - D 



As it follows from Theorem 0J spec{A,B) = {A: r(P D P cw ) = 0} = {A: r{h\) = 0}, 
where 

(44) h x (x) = (A + A 9 Bx) A (-X + B^Ax). 

The function h\ can be represented as infimum of max-plus linear maps so that the 
selection property f )37|) is satisfied. Namely, 

(45) h x (x) = AH^x, 

p 

where for i — 1 , . . . , n 



A - a ki + B k ., for 1 ^ k ^ n, a ki ^ -00 
-A - b k i + A k ., for 1 < k $ 

the brackets meaning that any listed choice can be taken 



(46) (tf?t , 

-A - b ki + Ak., for 1 < k ^ n, b ki ^ 
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The greatest eigenvalue of H\ equals the maximum cycle mean of H\. Using formula 
f )4T|) . we observe that r(H\) is a piecewise-affine function, meaning that it is composed of 
a finite number of affine pieces. More precisely, we have the following. 

Proposition 5. r(H^) is a finite piecewise-affine convex Lipschitz function of X. 

Proof. Using (j4Tj) we observe that r(H^) = — oo if and only if the associated digraph 
of H x is acyclic, which cannot happen when A and B and hence do not have — oo 
rows. 

If r(H^) is finite, then any finite cycle mean of H x can be written as {kX + a) /I, 
where / is the length of the cycle and k is an integer number with modulus not greater 
than /, hence this affine function is Lipschitz. The function r(H^) is pointwise maximum 
of a finite number of such affine functions, hence it is a convex Lipschitz piecewise-affine 
function. □ 

□ 

Definition 1 (Spectral Function). We define the spectral function of ( WZ\i by 

(47) s(X):=r(h x ) = r(P D P c{x) ). 

It follows from Theorem H] that s(X) ^ and that s(X) = if and only if A £ spec(v4, B). 
In general, s(X) is equal to the inverse minimal Chebyshev distance between Ax and 
A + Bx. 

By Proposition HJ 

(48) s(X) = Ar(H { x p) ). 

v 

As r(H ( x p) ) are piecewise-affine and Lipschitz, we conclude the following. 

Corollary 1. s(A) is a finite piecewise-affine Lipschitz function. 

Let us indicate yet another consequence of the fact that r(H ( x p) ) and s(A) are piecewise- 
affine. 

Corollary 2. 7/spec(A, B) is not empty, then it is a finite system of closed intervals and 
points. 

Note that this also follows, by means of projection, from a result by De Schutter and De 
Moor (30] that the solution set of a system of polynomial (in)equalities in the max-plus 
algebra is a (finite) union of polyhedra. The method of De Schutter and De Moor can 
also offer an alternative (computationally expensive) way to determine the spectrum and 
the generalized eigenvectors. 
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Conversely, it is shown in [4T] that any system of closed intervals and points in R can 
be represented as spectrum of (A, B). See also Subsect. 13.51 

3.2. Bounds on the spectrum of (A, B). Next we recall a bound on the spectrum 
obtained by Cuninghame-Green and Butkovic [151 121], extending it to the case when 
A = (ay) and B = (&y) may have infinite entries. Denote 

D(A,B)= V AijBi., 

(49) i: M- finite 

D(A,B) = - V Bi.fAi.. 

i: Bi. finite 

We assume that V = — oo and — V = +oo. 

Since Aj./Sj. = max{7 | A4. ^ 7 + B^.} is finite when A^. is finite and Sj. is not —00, 
we immediately see the following. 

Lemma 2. D_(A, B) (resp. D(A, B) ) is finite if and only if there exists ani e {1, . . . , m} 
such that Ai. is finite (resp. B^. is finite). 

When A and B have finite entries only, D_(A, B) and D(A, B) are just like the bounds 
of j2H Theorem 2.1]: 

D(A,B) = VA(0tf-6tf), 

(50) 

£>(A,fl) = AV(a ij -i ij ). 
* 3 

Note that £>( A B) and Z>(A, 5) defined by ([49]) take infinite values if A or 5 do not 
contain any finite rows. 

Proposition 6. If Ax ^ A + Bx (resp. Ax ^ A + Bx) has solution x > —00, then 
A > L>(A 5) fresp. A < D(A, B) ). 

Proof. If there exists i such that > A + 6^ for all j = 1, . . . , m, then Ax ^ A + Bx 
cannot have solutions. This condition is equivalent to A^jB^. > A plus the finiteness of 
A{.. Taking supremum of A^Bp over % such that A%. is finite yields D_(A,B). This shows 
that if Ax ^ A + Bx then A ^ D_(A, B). The remaining part follows analogously. □ 

□ 

The next result is an extension of [211 Theorem 2.1]. 
Corollary 3. spec(A,B) C [D(A,S),D(A,S)]. 
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We use identity (Tl3l) to give a more precise bound. It will be assumed that A and B 
do not have — oo columns. Note that this condition is more restrictive than that A and B 
do not have common — oo columns, and it cannot be assumed without loss of generality. 

Theorem 5. Suppose that A = (oy),5 = (6^) G Kmax™ ^° nc ^ ^o^e — oo columns. Then 

(51) spec(A, 5) C [-r(A s 5), r(E tt A)] C [D(A, B),D(A, B)\. 

Proof. Let Ax = XBx, then we also have 

Ax ^ X + Bx -X + x ^ A^Bx, 

(52) „ 

X + Bx^AxoX + x^ B*Ax. 

As A and B do not have — oo columns so that A^Bx and B^Ax do not have +oo entries, 
we can use 1115]) to obtain from (J32J that A G [-r(A tt S), r(M4)]. For A = r{B*A) we can 
find y ^ — oo such that A + y ^ B^Ay and hence A + -By ^ Ay. Using Proposition [6] we 
obtain A ^ D(A, B). The remaining inequality A ^ Z2(A, B) can be obtained analogously. 
□ 

□ 

By comparison with the finer bounds — r(A^B) and r(B^A), the interest of the bounds 
of Butkovic and Cuninghame-Green, D_(A, B) and D(A, B), lies in their explicit character. 
However, these bounds become infinite when the matrices A and B do not have any finite 
rows. We next give different explicit bounds, which turn out to be finite as soon as A and 
B do not have any identically infinite columns. 

Proposition 7. We have 

spec(A,B)C |J [-(AtB0)i, (B»iO),] , 

and so 

spec(A,B) C [- V(A*B0)i, V(B»A0)i], 

i i 

where zs £/ie n-vector of all 's. 

Proof. Consider x := and /i := Vj[/iA(0)]i, so that h\(x) ^ fi + x. Then, the non- 
linear Collatz-Wielandt formula ( fl4l) implies that r(h\) ^ If A G spec(A, 5), we have 
^ t{H\), and so, there exists at least one index % G {1, . . . , n} such that 

< [h x (0)]i = (A + (A«S0)i) A (-A + (fl<Uo) 4 ) . 

It follows that A ^ (B*A0)< and A ^ -(A*B0)i. □ 

□ 
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Figure 1 . Spectral function of f )53|) 
Remark 1. It follows readily from the Collatz-Wielandt property ( I14p that 

[-r(A«5),r(E tt A)] C [- V(A tt B0)i, V(5 tt A0)i] 

i i 

Example 1. We next give an example, to compare the bounds of Corollary [3], Theorem [5] 
and Proposition [7J Consider the following finite matrices of dimension 3x4: 

(-2 3 -3 -3\ /-4 5 -3 3 \ 



(53) 



.4 



B 



2 





2 



-1 
-3 



-v 



-4 1 2-2 
\5 -1 5 -1/ 

From the graph of spectral function, Figure [U it follows that the only eigenvalue is —2 since 
s(— 2) = and s(A) < for any A 7^ —2. The interval [— r(A^B), r(B^A)] is in this case 
[-2, 0.5]. Bounds ([50]) of [HI Theorem 2.1] yield the interval [D(A, B),D(A, B)] = [-3, 2], 
which is less precise. Proposition [7J yields the union of intervals [3, 0] = 0, [—2, —2], [3, 3] 
and [—3, —2], thus [—3, —2] U {3}. Note that these intervals are incomparable both with 
[-r{A*B),r(BtA)\ and \D(A, B), D(A, B)\ = [-3,2]. 

We remark that the intervals [-\Z i (A i B0) il \/ i (B i A0) i \ and [D(A, B),D(A, B)] are 
also in general incomparable. Also, Subsect. 13.51 will provide an example where the 
bounds [— r(A^B), r(B^A)] are exact. 
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Example 2. Let us now illustrate the discrete event systems interpretation of the spectral 
problem of the previous example. For readability, we replace the matrices by 



(54) 



A 



f-2 

— oo 

V 5 



3 
1 

-oo 



-oo 
2 
5 



-oo\ 
-oo 

- 1 / 



B 



( — oo 
2 



5 

-oo 
2 



-3 

-oo 
-oo 



-oo\ 
4 

-oo J 



This pair of matrices can be shown to have the same spectral function (Figure [TJ as the 
previous one, and the same bounds [— r(A^B), r(B^A)}. Consider now the two discrete 
event systems 

y = Ax, z = Bx . 

Here, Xi is interpreted as the starting time of a task i, and yi and Zi are interpreted as 
output time. This is illustrated in Figure [2j For instance, the constraint y\ = max(— 2 + 
x i,3 + X2) in y = Ax expresses that the first output is released at the earliest, given 
that it must wait 3 time units after the second input becomes available, and can not be 
released more than 2 time units before the first input becomes available. We are looking 
for a common input x such that the time separation between events is the same for both 
outputs, so that 



Vi ~ Vj = Zi 



This can be solved by finding an eigenvector x, so that Ax = A + Bx. By inspection of 
the spectral function in Figure [fl we see that A must be equal to —2. Then, computing x 
reduces to solving a mean payoff game (see the discussion in section 13.41 below for more 
algorithmic background). In this special example, x can be determined very simply by 
running the power type algorithm (like the alternating method of [20] ) 



x 



(0) 



(0,0,0,0) T , x (fe+1) = h_ 2 (x ik) ) 



where 



h_ 2 (x) := (-2 + A*Bx) A (2 + B^A, 



x 



until the sequence x k converges. Actually, 



x 



(2) 



~( 3 ) = (_5n,-5, 



and it can be checked that 



^ 2 ) = -2 + J Bx( 2 ) = (3,l,0) T 




FIGURE 2. Finding a common input making the outputs of two discrete 
event systems indistinguishable, modulo a constant 

3.3. Asymptotics of the spectral function. If A and B do not have — oo columns, 
the functions X + A^B and — X + B^A are represented as infima of all max-linear mappings 
and, respectively, such that 

{K^ ] )i. = X - a ki + B k ., 1 ^ k < n, a ki ^ -oo, 

( 55 ) i s 

{m£% = -X - + A fc ., 1 ^ k ^ n, b kl ^ -oo. 

This representation satisfies the selection property. 

Matrices and m[ s) are both instances of H^' which represent h\. We will need 
the following observation on r(Hjjf') 

Lemma 3. Denote k := min(2m, n). The spectral radii r(H^) can be expressed as 
Xs/l + a, where ^ \s\ ^ / ^ k, and a ^ A (A, B), where 

(56) A(A,S):= V (a i3 -b ik )V \J (k 3 -a ik ). 

i,j,k: ciijjt-oo, b ik ^-oo i,j,k: 6y^-oo, a^^-oo 

Moreover it is only possible that s = I — 2t for t = 0, . . . , I. 

Proof. According to ff¥Tl) and (HEl) . r(H^) is a cycle mean of the form 

(57) (*&, + ••• + &)/* 

where we use the notation t\ 3 - := —a k i + X + b k j and = —A — b k % + a k j for i, j = 1, . . . , n 
and k = l,...,m. Actually £f is the entry of but here we also need the 

intermediate index k. Note that it is determined by i. 

In (ISTj) . only even numbers of ±A can be cancelled, hence it can be expressed as Xs/l + a 
where ^ \s\ ^ I with s = I — 2t for £ = 0, . . . , I. The cycle . . . , ij) is elementary, 
hence Z ^ n. We also obtain a ^ A (A, 5) since the arithmetic mean does not exceed 
maximum. 

It remains to show that I ^ 2m. Indeed if I > 2m then there is an upper index which 
appears at least twice in (JoTj) . Assume w.l.o.g. that this is k±. Then the sum in (|57p takes 
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one of the following forms: 

— + [A + bk li2 + . . . — a klir ] + A + bk r i r+1 + ■ ■ ■ , 

(58) 

— A — fefciii + [flfcii 2 + ... — A — 6fcii r ] + a fcii r+ i + • • • 

Assume w.l.o.g. that we have the first one. Then we can split it into the following two 
cycles, as indicated by the square bracket in the first line of (1581) : 

±kx _|_ ±ki _|_ Jtr- 1 

V«2 ' «2«3 ' ' ' ' V-lV' 

(59) 

-fel _|_ .fer + 1 _|_ .fcz 

H«r + 1 «r+l«r+2 ' * ' Ml' 

Each of these forms is a weight of a cycle in Indeed, (1581) (the first expression) 

indicates that k\ is chosen by i r in H^y so that any element t*^ for j = 1, . . . , n is an 
entry of H^. All other elements in (159]) are also entries of if^. 

The arithmetic mean for both of the cycles in (15 9 p has to be equal to (15 7p . if this is 
indeed r(H^). This shows I ^ 2m. □ 

□ 

We also define 

C(A,B):= \/ faj -*>*), 

i.j.k: aaj^—oo, bu-j^—oo 

(60) 

C(A,B):= f\ (Oik-bij). 

i,j,k: a ik ^-co, by^-oo 

We now study the asymptotics of s(A), both in general case and in some special cases. 

Theorem 6. Suppose that A,B£ IR^^T and denote k := min(2m, n). 

1. There exist k\, l\, k 2 , l 2 such that ^ l\ ^ k, k\ = li — 2t\ where ^ t\ ^ [A/^J? 
^ l 2 ^ k, k 2 = l 2 — 2t 2 where ^ t 2 ^ L^/2J, and ai,a 2 ei sttc/i £/iai 

s(\) = \k 1 /h + a 1 , if\^-2K 2 A(A,B), 

(61) 

s(A) = -Xk 2 /l 2 + a 2 . if 2k A(A, B). 

2. Suppose that A and B do not have —00 columns. Then there exist a\ ^ r(A^B) 
and a 2 ^ r(B^A) such that 

s(\) = \ + a 1 , ifX^-2KA(A,B), 

(62) 

s (A) = -A + a 2 , if X^2kA(A,B). 

3. Suppose that A and B are real. Then 

s(X) = X + r(A i B), ifX^C(A,B), 

( 63 ) ^ « - 

s(A) = -A + r(BU), ifX^C(A,B). 
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Proof. 1: For the proof of this part, we observe that for each A, the function s(A) is the 
maximum cycle mean of a representing matrix so that it equals Xk/l + a where 

^ / ^ k, k = I — 2t where ^ t ^ I. For any two such terms, difference be- 
tween coefficients kjl is not less than 1/k 2 , and the difference between the offsets does 
not exceed 2A(A, B), which yields that all intersection points must be in the interval 
[-2k 2 A(A 1 B),2k 2 A(A,B)}. Thus s(A) is just one affine piece for A ^ -2k?A(A,B) and 
for A ^ 2k, 2 A(A,B). As s(A) ^ for all A, the left asymptotic slope is nonnegative, and 
the right asymptotic slope is non-positive. 

2: When A does not have — oo columns, some of the matrices are of the form 
and their maximum cycle mean is A + a. Taking minimum over all r(H^) of the form 
A + a yields an offset at ^ r(A*B). The cycle mean A + a>\ will dominate at small A, 
and the smallest intersection point may occur with a term A(k — + a[. Indeed, the 
difference between coefficients is precisely the smallest possible and the difference 
\a\ — a[\ may be up to 2A(A,B). This yields the bound — 2kA(A, B). An analogous 
argument follows when A is large and B does not have — oo columns. 
3: When A and B are real and A < C_(A,B), all coefficients in the min-max function 
A + A^B are real negative, and all coefficients in the min-max function —A + B^A are real 
positive. This implies that s(A) is equal to the minimum over r(K^), which is equal to 
A + r(A^B). An analogous argument follows when A > C(A,B). □ 



In Proposition [TU1 we will show by an explicit construction that any slope kjl can be 
realized as asymptotics of a spectral function. 

We next observe that the asymptotics of s(A) can be read off from the spectral function 
s°(A), which we introduce below. For arbitrary C = (c^) G M™^ x n define 



Let s°(X) be the spectral function of the eigenproblem A°x = A + B°x. 

Proposition 8. Suppose that A, B E IRmax™ an d ^at \k\jl\ where k\,l\ ^ (resp. 
—\k2ll2 where k2,h ^ 0) is the left (resp. the right) asymptotic slope of s(X). Then 



□ 



(64) 




(65) 
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Proof. Observe that the representing matrices of 



(66) 



h° x := (A + (A°)»5°a;) A (-A + {B°fA°x) 



are in one-to-one correspondence with the representing matrices of h\. The finite 
entries equal to ±A, they are in the same places and with the same sign of A as in 

Hyf 1 . Hence the cycle means in have the same slopes as the corresponding cycle 

means m H%\ but with zero offsets. When s(A) = r{h\) is computed by (j4"5l) . the 
asymptotics at large and small A is determined by the slopes only and yields the same 
expression as for s°(A) = r(h° x ). □ 



3.4. Mean-payoff game oracles and reconstruction problems. Here we consider 
the problem of identifying all afline pieces that constitute the spectral function and com- 
puting the whole spectrum of (A, B) in the case when A and B have integer entries. 

The result will be formulated in terms of calls to a mean-payoff game oracle (computing 
the value of a mean payoff game). Let us briefly describe what the mean-payoff games are 
and how they are related to our problem. For more precise information the reader may 
consult Akian et al. [5] and Dhingra, Gaubert [23], as well as Bjorklund, Vorobyov [TTJ 
and Zwick, Paterson |44j . 




sum deterministic mean-payoff game, which also corresponds to the system Ax ^ Bx. A 
schematic example of such a game is given in Figure El left. Two players, named Max 
and Min, move a pawn on a bipartite digraph, whose nodes belong either to Max (□) 
or to Min (O)- I n the beginning of the game, the pawn is at a node j of Min, and she 
has to move it to a node i of Max, paying to him — ay (some real number). Then Max 
has to choose a node k of Min. While moving the pawn there, he receives from her. 
The game proceeds infinitely long, and the aim of Max (resp. Min) is to maximize (resp. 
minimize) the average payment per turn (meaning a pair of consecutive moves of Min 
and Max). It turns out that the game has a value, which depends on the starting node 
of Min. Moreover r(A^B) equals the greatest value over all starting nodes (i.e., all nodes 
of Min). 

The two-sided eigenproblem Ax = A + Bx can be represented as 



□ 




It can be observed that the min-max function A^B is also a dynamic operator of a zero 



(67) 
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Figure 3. General mean-payoff game (left) and mean-payoff game corre- 
sponding to Ax = A + Bx (right) 

This is equivalent to x ^ h\(x) where h\(x) := (A + A^Bx) A (—A + B^Ax) as above. 
Hence the problem Ax = A + Bx corresponds to a parametric mean-payoff game of special 
kind, with 2m nodes of Max and n nodes of Min, whose scheme is displayed on Figure |3j 
right, where individual nodes of the players are merged in three large groups. 

Denoting by MPG (m,n,M) the worst-case execution time of any mean-payoff oracle 
computing r(A^B), where A,Bg M™ x n have — oo entries and integer entries with the 
greatest absolute value M, we immediately obtain that for the same A and B we can find 
s(0) = r(h) by calling that oracle, in no more than MPG(2m, n, M) operations. 

The implementation of a mean-payoff oracle can rely on the policy iteration algorithm 
of [El [23] , as well as the subexponential algorithm of [H] or the value iteration of jB] . 
Zwick and Paterson jH] showed that MPG(m, n, M) is pseudo-polynomial. We use this 
result below to demonstrate that the graph of spectral function s(A) can be reconstructed 
in pseudo-polynomial time. 

Theorem 7. Let A,B G R^ax™ have only — oo entries and integer entries with absolute 
value bounded by M. Denote k := min(2m,n). 

1. All affine pieces that constitute the function s(A) and hence the spectrum of (A, B) 
can be identified in no more than A(A, B)0(k 6 ) calls to the mean-payoff game 
oracle, whose worst-case complexity is MPG(2m, n, k, 2 (M + 4Mk 2 )). In particular, 
the reconstruction can be done in pseudo-polynomial time. 

2. When A and B have no — oo columns, the number of calls needed to reconstruct 
the function s(A) can be decreased to A(A, B)0(k 5 ), where each call takes no more 
than MPG(2m, n, k 2 (M + 4Mk) operations. When A and B are real, the number 
of calls is decreased to (C(A, B) — C_(A, B))0(k a ), and the complexity of each call 
to MPG(2m, n, 3Mk 2 ) operations. 
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Proof. In all cases we have a finite interval L of reconstruction, determined by the asymp- 
totics of s(A). Using Theorem [6], we obtain that in case 1 this is 

(68) L := [-2k 2 A(A, B), 2k 2 A(A, B)] C [-4k 2 M, 4k 2 M], 
In case 2, this is 

(69) L := [-2«A(A, 5), 2kA(A, 5)] C [-4«M, 4kM] 
when A and -B do not have — oo columns, or 

(70) L := [C(A,B),C(A,B)} C [-2M, 2M] 

when A and -B do not have — oo entries. 

We first compute the asymptotic slopes of s(A) outside L. By Proposition [HJ we can 
do this by computing s°(±l) in just two calls to the oracle which computes it in no more 
than MPG(2m, n, 1) operations. Then the goal is to reconstruct all affine pieces which 
constitute s(A) in the interval L. 

The affine pieces of s(A) correspond to the maximal cycle means in the matrices from 
the representation of h\(x). The points where such affine pieces may intersect are given 
by 

, ai + k\X _ a 2 + k 2 \ 

where all parameters are integers and 1 ^ \ki\, l^li^ii^ ^ k by Lemma [31 This implies 

(72) A - aiU2 - ° 2ni 

k 2 n 1 - kin 2 

The denominators of these points range from — n 2 to k 2 , hence their number is \L\0(k 4 ) 
where \L\ is the length of the reconstruction interval L. We reconstruct the whole spectral 
function by calculating s(A) at these points, since there is only one affine piece of s(A) 
between them. 

Using ( )68|) . ( )69|) and (!70|) we obtain that the absolute value of the entries of A and 
A + B at each call does not exceed M + 4k 2 M in case 1, and M + 4kM or M + 2M in case 
2. Multiplying the entries of A and A + B by the denominator of A which does not exceed 
k 2 , we obtain a problem with integer costs, where all maximum cycle means r(H^) get 
multiplied by that denominator, and hence s(A) gets multiplied by that denominator as 
well. Thus we can solve this mean-payoff game instead of the initial one. In case 1, the new 
integer problem can be resolved by the mean-payoff oracle in MPG(2m, n, k, 2 (M + 4k 2 M)) 
operations. In case 2, it takes no more than MPG(2m,n, k 2 (M + 4kM)) operations when 
A and B do not have — oo columns, and no more than MPG(2m, n, 3Mk 2 ) operations 
when A and B do not have — oo entries. The proof is complete. □ 
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□ 

Since spec(A,B) is the zero set of s(A), we can identify spec(A,B) by reconstructing 
s(A) in the intervals given by Proposition [7] or more generally, Theorem [6j However, the 
task of reconstructing spectrum of (A, B) as zero-level set is even more simple, by the 
following arguments. 

Theorem 8. Let A, B e H^ax" have only integer or — oo entries. 

1. In general, the identification of spec(A, B) requires no more than MO(k 3 ) calls to 
the mean-payoff' game oracle, whose worst-case complexity is MPG(2m, n, 2k(M + 
2Mk)). In particular, spec(A, B) can be identified in pseudo-polynomial time. 

2. If A and B have no — oo columns, then the number of calls to the oracle needed 
to identify spec(A, B) does not exceed (V i ( J B tt A0) i + V i (A tt J B0) i )O(K 2 ), and the 
complexity of the oracle does not exceed MPG(2m, n, 6Mk) operations. 

Proof. We have to reconstruct the zero-level set of s(A), within a finite interval L of 
reconstruction. In case 1, we notice that the intersection of s(A with zero level can occur 
only at points with absolute value not exceeding 2Mn (since s(A) consists of affine pieces 
(a + kX)/l where \a\ ^ 2Mk). Hence in case 1 

(73) L := [—2Mk, 2Mk]. 
In case 2 we use the bounds of Proposition [7J 

(74) L := [- V(A tt S0)i, V(S s A0) i ] C [-2M, 2M] 

i i 

when A and B do not have — oo columns. In case 1, we also need to check the asymptotics 
of s(A) outside the interval, for which we check s°(±l) = (i.e., s°(±l) ^ which takes 
no more than MPG(2m, n, 1) operations). 

The absolute value of entries of A and A + B does not exceed M + 2Mk in case 1 and 
M + 2M in case 2. We have to check s(A) = (i.e., s(A) ^ 0) at all possible intersections 
of affine pieces constituting s(A) with zero, i.e., at the points A = a/k within L, such that 
a and k are integers and k ^ k. We also may have to check s(A) = for one intermediate 
point between each pair of neighbouring points Ai and A 2 such that s(Ai) = s(A 2 ) = 0. If 
it holds then s(A) = holds for the whole interval, and if it does not then it holds only 
at the ends. Note that such an intermediate point for ax/k\ and a 2 //c2 can be chosen as 
(ai + a 2 )/(fci + k 2 ) thus leading to k ^ 2k. 

Multiplying all the entries by k yields a mean-payoff game with integer costs, for which 
we check whether the value is nonnegative. This takes no more than MPG(2m, n, 2k(M + 
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2Mk)) in Case 1 and MPG(2m,n,2ft x 3M) in Case 2, with the number of calls not 
exceeding \L\0(n 2 ). □ 

□ 

Note that this theorem uses the oracles checking s(A) is 0, not requiring to compute 
the exact value. 

We can also formulate a certificate that A is an end (left or right) of a spectral interval. 

Proposition 9. Supose that s(A*) ^ 0. Then A* is the left (resp., the right) end of an 
interval of spec(A, B) if and only if there exists a representing matrix where the 
weights of all cycles are nonpositive, and the slopes of all cycles with zero weight are 
strictly positive (resp., negative). 

Proof. Condition s(A*) ^ assures that A* G spec(v4, B). We will consider the left 
end case, for the right end the argument is similar. Recall that s(A*) := r(h\* admits 
an inf-representation with selection property. Since there is only finite number of 
representing matrices, there exists such that s(A) = r(H^ for A G [A* — e, A*] (for 
the right end, we would consider A G [A*, A* + e]). Then A* is the left end of a spectral 
interval if and only if r(H^) = but r(H^) < for A G [A* — e, A*). After applying the 
definition of f l46|) and the maximum cycle mean formula for r(H^) (HT|) . the claim 
follows. □ 

□ 

Observe that the condition in Proposition |9] can be verified in polynomial time for a 
given Namely, it suffices to compute the maximum cycle mean, identify the critical 

subgraph consisting of all cycles with zero cycle mean, and solve the maximum cycle mean 
problem for that subgraph, with the edges weighted by 1 or —1 according to the choice 
of A or -A in (|46l) . 

The reconstruction of spectral function has been implemented in MATLAB, also to 
generate Figures [TJ and HI 

3.5. Examples of analytic computation. In this section we consider two particular 
situations when the spectral function can be constructed analytically. The first example 
shows that any asymptotics k/l, where / = 1, . . . , m and k = I — It for t = 1, . . . , I, 
can be realized. The second example is taken from [H] , and it shows that any system 
of intervals and points on the real line can be represented as spectrum of a max-plus 
two-sided eigenproblem. 
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Asymptotic slopes. In our first example we consider pairs of matrices A m ' 1 £ 
MSSTN-B"*' 1 e M™ a x x m with entries in {0, -00} , where ^ / ^ [m\. An intuitive idea 
is to make some "exchange" between the max-plus identity matrix and some cyclic per- 
mutation matrix. For instance 



(75) 













/o 




5 6 








\ 



0/ 



where the dots denote —00 entries. 



are defined as matrices with {0, — 00} entries such that a™' = 



m,l 

m, or 2 = j + 1 where 21 < i ^ m, or i = j = 2k where 1 ^ k ^ I, or 



m,( 



2&, where 1 ^ < /, and a^' 

1 m,l\ 



-00 otherwise. 



Formally, A 
for 2 = 1 and j = 
i = 2k + 1 and j 

Similarly, B m ' 1 = (b^' 1 ) are defined as matrices with entries in {0, —00} such that 
b™' = for i = j where 21 < i ^ m, or i = j = 2k — 1 where 1 ^ k ^ I, or i = 2k and 
j = 2k — 1, where 1 ^ k ^ I, and = —00 otherwise. 

Proposition 10. The spectral function associated with A m, \B m ' 1 consists of two linear 
pieces: s(A) = A(m — 21) /m for A ^ and s(A) = — A(m — 21) /m for A ^ 0. 



Proof. Let us introduce yet another matrix C' m ' (A) 



m,Z 



(A)) G 



omxm 
v max " 



Informally, it 



is a sum of a {0, —00} permutation (circulant) matrix and its inverse, weighted by ±A. 
This pattern corresponds to the above mentioned "exchange" in the construction of A m ' 1 
and B m ' 1 . In particular, f lTB]) corresponds to 



(■ 

A 



(76) 



C 



6.2 



-A 



-A 



-A 



-A 



Defining formally, c[ 



1 11, 1 

m 



-A, c 



m,l 



A, and 



(77) 



ml 



sign(z, j)A, if 1 ^ i,j ^ m and |j — i\ 
-00, otherwise, 
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where 

(78) sign(i,j) 



1, j - 1 = % > 21 or j ± 1 = % = 2k < 21, 
-1, z - 1 = j ^ 2/ or % ± 1 = j = 2k < 2/. 



Observe that the pairs (i, j) and (j, i) for j = i + 1 (and also (l,m) and (m, 1)) have the 
opposite sign. 

It can be shown that each representing max-plus matrix of the min-max function 
(79) K'\ x ) = ( A + (A m ' l fB m ' l x) A (-A + (B m ' l fA m ' l x) 

is choosing one of the two entries in each row of C m > l (\). The matrices can be classified 
according to this choice as follows (see (1761) for example): 

1. Choose (m, 1), and (i, i + 1) for i — 1, . . . , m — 1; 

2. Choose (1, m), and (i, i — 1) for i = 2, . . . , m; 

3. Choose both (m, 1) and (1, m), or both (i — 1, i) and (i, z — 1) for some % — 2, . . . , n. 
The first two strategies give just one matrix each, with the (maximum) cycle means 
A(m — 21) /m and — A(m — 21) /m. The rest of the representing matrices are described by 
3., and it follows that their maximum cycle means are always greater than or equal to 0. 
Hence s(A) = A(m - 21) /m A -A(m - 21) /m. □ 



□ 



The spectrum of two-sided eigenproblem. Now we consider an example of 
Let us define A e R^, B G R 2 ™: 



2a,- 26, 2c,- 




.4 

10) 

B 



where ^ q < a i+1 for i = 1, . . . , t — 1, where 6, := 2i±s. The following result describes 
spec(^4, B). 

Theorem 9 (|3Tj). With A, B defined by (IHUl) . 

t 

(81) spec(AS) = |jKcj]- 

To calculate s(A), which is a more general task, one can study the representing matrices 
like in the previous example. Another way is to guess, for each A, a finite eigenvector of 
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1.5 




Figure 4. The spectral function of A and B in 



PdPc(\) and then s(A) is the corresponding eigenvalue. By this method we obtained that: 



(82) 



s{\) = { 



A — ai, if A ^ d, 

0, if a k ^ A ^ c k , k — 1, . . . , t, 

max(c fc - A, A - a k+1 ), if c k ^ A ^ a fc+ i, fc = 1, — , t — 1, 

Q — A, if A ^ q. 



More precisely, it can be shown that the following vectors are eigenvectors of PdPc(X)'- 



(83) y } 



(0 ai ai), if A ^ ai, 

(0 X + b k -a k A + 6 fc -a fc ), if a k < A < k = 1, . . . , t, 

(0 c fc c fc ), 

(0 A A), 

k (0 ct q) t , 



if 6 fe < A < c fe , k = 1, . . . ,t, 



if c fc < A < a fc+ i, fc = 1, . . . , t - 1, 
if A ^ q, 



with the eigenvalues expressed by ( 1821 . 

We can also conclude that in this case —r(A^B) = a x and r(B*A) = c t . Indeed, by 
( 1521 . s(A) = A — ai for A ^ ai and s(A) = c t — A for A ^ q. Comparing this with the 
result of Theorem [6], part 3, we get the claim. 

As a\ and c t are eigenvalues, the last result shows that the bounds given in Theorem [5] 
cannot be improved in general. 
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For example, take t = 3, [a\, Ci] = [1, 2], [a 2 , c 2 ] = [2.2, 2.4] and [a^, C3] = [3, 3]. 
Then 

_ /l 1.5 2 2.2 2.3 2.4 

" \ 2 3 4 4.4 4.6 4.8 

(84) ) 

5 _/ooo 

\1 2 1.5 2.2 2.4 2.3 

The spectral function is shown on Figure HI Note that this is the least Lipschitz function 
with a given zero-level set. The same observation holds for the general case (|82|) . 

4. Conclusions 

We have developed a new approach to the two-sided eigenproblem A®x = A Cg> £> Cg> x in 
max-plus linear algebra, based on parametric min-max functions. This yields a reduction 
to mean-payoff games problems, for which a number of algorithms have already been 
developed. We introduced the concept of spectral function s(A), defined as the greatest 
eigenvalue of the associated parametric min-max function (or the greatest value of the 
associated mean-payoff game). We showed that s(A) has a natural geometric sense being 
equal to the inverse of the least Chebyshev distance between A <g> x and B <g> x. The 
spectrum of (A, B) can be regarded as the zero-level set of the spectral function, which 
is a 1-Lipschitz function consisting of a finite number of affine pieces. These pieces can 
be reconstructed in pseudopolynomial time, hence the spectrum of (A, B) can be also 
effectively identified. 

A similar approach can be used in max-plus linear programming [27J. Spectral functions 
of a different type are used in the decision procedure associated with the tropical Farkas 
lemma in [Sj, allowing one to check whether a max-plus inequality can be logically deduced 
from other max-plus inequalities. The present approach can be generalized to the case 
when the entries of A and B are general piecewise-affine functions of A |43j . but the case 
of many parameters would be even more interesting. Such development could lead to 
practical applications in scheduling and design of asynchronous circuits. Also note that 
the parametric tropical systems are equivalent to parametric mean-payoff games, directing 
to useful stochastic and infinite-dimensional generalizations. 
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